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Abstract. Tunnelling processes through black hole horizons have recently been investigated in the 
framework of WKB theory, discovering an interesting interplay with the Hawking radiation. In this 
paper, we instead adopt the point of view proper of QFT in curved spacetime, namely, we use a suitable 
scaling limit towards a Killing horizon to obtain the leading order of the correlation function relevant 
for the tunnelling. The computation is done for certain large class of reference quantum states for scalar 
fields, including Hadamard states. In the limit of sharp localization either on the external side or on 
opposite sides of the horizon, the quantum correlation functions appear to have thermal nature. In both 
cases the characteristic temperature is referred to the surface gravity associated with the Killing field and 
thus connected with the Hawking one. Our approach is valid for every stationary charged rotating non- 
extremal black hole. However, since the computation is completely local, it covers the case of a Killing 
horizon which just temporarily exists in some finite region, too. These results provide strong support 
to the idea that the Hawking radiation, which is detected at future null infinity and needs some global 
structures to be defined, is actually related to a local phenomenon taking place even for local geometric 
structures (local Killing horizons), existing just for a while. 



1 Introduction 

As is known, the Hawking radiation [H a75| is detected at future null infinity of a spacetime 
containing collapsing matter giving rise to a black hole. At least in the case of spherical sym- 
metry, the existence and the features of that radiation are quite independent from the details 
of the collapse. However, the type of the short-distance behaviour of the two-point function 
of the reference state, employed to describe the modes of the radiation, plays a relevant role 
[FH90]. In recent years, attention has been focused on local properties of models where the 
Hawking radiation is manifest. Here local means in a neighborhood of a point on the event 
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horizon |PW00j lA"NVZ05j lA"PgQ8l IAPGS08] . In this second approach the radiation appears to 
be related with some thermal effects associated to some tunnelling process through the horizon. 
In particular, the tunneling probability, computed in the framework of semiclassical WKB ap- 
proach, has the characteristic thermal form e~ E ^ Ta (the Boltzmann constant being re-defined 
as k = 1) where Tjj is the Hawking temperature and E is the energy of the particle crossing the 
horizon. More precisely, that exponential thermal factor arises when taking a limit towards the 
horizon for an endpoint of the path of the classical particle. 

This approach is interesting because it deals with local aspects only and, in this sense, it seems 
to be more general than the standard one. In fact, it may be applied to pictures where a certain 
geometric structure, interpreted as the horizon, exists "just for a while", without extending 
into a true global structure up to the future null infinity where, traditionally, the Hawking 
radiation is detected. In [DNVZZ07, HDVNZ09], even the case of a spherically symmetric black 
hole in formation was analysed, where no proper horizon structure exists, being replaced by a 
dynamical horizon. Other interesting results, also considering the backreaction, can be found in 
[KM071IMV05] . 

Within these new remarkable approaches it is however difficult to understand how the found 
properties are independent form the state of the quantum system. This is essentially due to 
the fact that they are discussed at the quantum-mechanical level rather than the quantum-field- 
theory level. Indeed, in [P W001 IANVZ05] it is assumed that there is some preferred notion 
of quantum particle whose wavefunction satisfies the Klein-Gordon equation. That equation 
is actually treated as a Schrodinger equation when dealing with transition probabilities within 
the WKB framework. However, in curved spacetime there is no natural definition of quan- 
tum particle, unless adopting the quantum-field-theory framework, fixing a preferred quasifree 
reference state and building up the associated Fock-Hilbert space. It does not seem that this 
relevant issue is tackled in the mentioned literature. Furthermore, the procedure exploited in 
[ANVZ05| IDNVZZ07, IHDVNZ09] needs a Feynman prescription to make harmless a divergence 
that pops up when performing the above-mentioned limit towards the horizon. As a matter 
of fact, that procedure turns the real-axis computation into a complex-plane computation and 
the very imaginary part of the WKB amplitude, arising that way, leads to the wanted factor 
e -E/T H _ reason f or appearance of an imaginary part in the semi-classical action has 
been clarified in [APS08, APGS08J, through a very careful analysis of the integration paths and 
the adopted coordinate system, without requiring any regularization procedure and remaining 
in the general WKB framework, while adopting the path integral viewpoint. In any case, that 
clarification does not solve the problem of the absence of a precise notion of a quantum particle 
in curved spacetime, necessary to exploit the WKB formalism at quantum-mechanical level. It 
seems however plausible that the appearance of an imaginary part in the semi-classical action, 
formally equivalent to the Feynman regularization procedure, is nothing but a remnant of the 
choice of a preferred reference state at quantum-field-theory level, whose two-point function has 
a short-distance divergence close to that of Minkowski vacuum. 

While sticking to the local aspects associated with states showing the Hawking radiation, dif- 
ferently from the references quoted above, in this paper we shall deal with a definite framework 
at the quantum-field-theory level. More precisely, we shall focus on the two-point correlation 
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function w(<&(cc)<&(?/)) of a quantum field $ settled in a (not necessarily quasifree) state u whose 
short-distance divergence is, essentially, of Hadamard type, thus generalising the short-distance 
behaviour of Minkowski vacuum. This is one of the hypotheses exploited in FIHK). 1KW91] . It 
actually encompasses a huge class of states, those that are supposed to have a clear physical 
meaning }Wa94j , especially in relation with the problem of the renormalisation of the stress- 
energy tensor and the computation of the quantum backreaction on the metric. 
The two-point function uj(<&(x)<&(y)) corresponds, up to normalization, to a probability ampli- 
tude. In this sense, it measures the tunnelling probability through the horizon when x and y are 
kept at the opposite sides of the horizon. In particular, if the state u is quasifree, that probabil- 
ity amplitude can be interpreted as a quantum-mechanical probability amplitude, as wished in 
the above-mentioned literature, just referring to the natural notion of particle associated with 
the state u. 

From the geometric viewpoint, we shall assume to work in a sufficiently small neighbourhood 
of a local Killing horizon structure "K, also supposing that the surface gravity k is nonvanishing 
and constant along the horizon. We stress that the structure could be either part of the future 
Killing event horizon of a stationary black hole in the full Kerr-Newman family also obtained by 
matter collapse, or it could be completely local and ceasing to exist in the future of in view of 
the general dynamics of the matter and the fields in the considered spacetime. The requirement 
that the surface gravity is constant on the local horizon means that, at least locally, a thermal 
equilibrium has been reached, since a constant surface gravity corresponds to the validity of the 
zero-law of black-hole thermodynamics. 

The existence of a timelike Killing vector K defining "K provides the preferred notions of 
time and energy we intend to consider. Notice that in [DNVZZ07, HD VNZ09] . dealing with 
spherically-symmetric black holes in formation, the notion energy was referred to the so called 
Kodama-Hayward vector field that, in those backgrounds, generalises the notion of Killing field. 

Exploiting general technical achievements about Killing horizons established in [KW91] and 
[RW92], we shall prove that, independently from the choice of the quantum state in above- 
mentioned class, when the supports of the test functions centred on the two arguments x, y of 
uj(&(x)&(y)) become closer and closer to the horizon, the two-point function acquires a thermal 
spectrum with respect to the notion of time and energy associated with the Killing field. More 
precisely, if both arguments stay on the same side of the horizon, the Fourier transform of the 
two-point function presents the very Bose-Einstein form driven by the Hawking temperature 
(see however the remark at the end of the paper). Conversely, whenever the two arguments 
are kept at the opposite sides of the horizon, the resulting spectrum is different. It is however 
in agreement with the transition probability between two weakly coupled reservoirs which are 
in thermal equilibrium at Hawking temperature. Actually, as in the case of the Boltzmann 
distribution, its spectrum decays exponentially for high energies. In both cases, in order to 
catch the leading contribution to the two-point function, we shall exploit a suitable scaling limit 
procedure [HNS84I IBu96[ IBV95] towards the horizon. Operating in this way, the local thermal 
nature of the correlation functions becomes manifest as a state-independent feature when the 
states belong to the above-mentioned wide class. 

The paper is organized as follows. In the next section, recalling some technical results 
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established in [KW91J and [RW92J, we shall present the geometric hypotheses we shall use. We 
assume that the reader is familiar with the standard notions of differential geometry of spacetimes 
[Wa84]. In the subsequent section we shall compute the two-point function u)(<&(x)<&(y)) and 
its limit approaching the horizon. The last section will present a summary and some general 
remarks. 

2 Spacetime Geometry 

2.1. Local geometry. We start our discussion fixing the basic geometric setup that we employ 
in this paper. We henceforth consider a 4-dimensional (smooth) time-oriented spacetime (M,g). 
Furthermore, we assume the validity of the following local geometric properties, which are the 
same as in jRW92j. 

Definition 2.1. Let be an open subset of M, the local general geometric hypotheses 

hold in if a smooth vector field K exists thereon such that: 

(a) K is a Killing field for g in 0. 

(b) contains a connected 3-submanifold K, the local Killing horizon, that is invariant 
under the group of local isometries generated by K and K a K a = on K. 

(c) The orbits of K in are diffeomorphic to an open interval contained in M. and K admits 
a smooth 2-dimensional cross section which intersects each orbit of K exactly once. 

(d) The surface gravity - i.e. the function k : IK — > R such that, in view of (a) and (b) 
^J a {KbK b ) = —2KK a - turns out to be strictly positiv^ and constant on K. 

As we said above, the local Killing horizon K may represent a horizon which exists "just for 
a while", without extending into a true global structure which reaches the future null infinity. 
However, our hypotheses are, in particular, valid [R W92| in a neighbourhood of any point on a 
black hole future horizon, once that, after the collapse, the metric has settled down to its sta- 
tionary, not necessarily static, form of any non-extreme black hole in the charged Kerr-Newman 
family. In particular, those hypotheses and our results are valid for the Kerr black hole. There 
K is the Killing vector defining the natural notion of time in the external region of the black 
hole and K is part of the event horizon. 

With the hypotheses (a) and (b), the integral lines of K along K can be re-parametrized to 
segments of null geodesies and \7 a (Ki,K b ) = —2KK a holds on IK where the surface gravity, 
k : K — > R, is constant along each fixed geodesic [Wa8j4j. The requirement (d) is not as strong as 
it may seem at first glance. Indeed, it is possible to prove that, whenever a spacetime admitting 
a Killing horizon satisfies Einstein equations and the dominant energy condition is verified, the 
surface gravity must be constant on the horizon [Wa84]. (However, \Va8 II a result originally 
obtained by Carter states that - independently of any field equations - the surface gravity of 
a Killing horizon must be constant if (i) the horizon Killing field is static or (ii) there is an 
additional Killing field and the two Killing fields are 2-surface orthogonal.) At least in the case 

1 What actually matters is k =^ 0, since k > can always be obtained in that case by re-defining K — > —K. 
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k > 0, this is nothing but the zero-th law of black hole thermodynamics where ^ amounts to 
the Hawking temperature of the black hole. 

2.2. Killing and Bifurcate Killing horizons. We now focus on the relation of the previously 
introduced local geometric hypotheses and the more rigid case of a bifurcate Killing horizon. A 
Killing field K determines a bifurcate Killing horizon [Bo69| when it vanishes on a connected 
2-dimensional acausal space-like submanifold 23 C M, called the bifurcation surface, and K 
is light-like on the two K- invariant 3-dimensional null submanifolds Jf + , !K_ C M generated 
by the pairs of null geodesic orthogonally emanated from 23. In particular "K+ n !K_ = 23 and 
the null geodesies forming "K + U "K- are re-parametrised integral lines of K on U \ 23. 
By definition, on "K+ the field K in the future of 23 is outward pointing. The simplest example 
of a bifurcate Killing horizon is that realized by the Lorentz boost K in Minkowski spacetime. 
Other, more interesting, cases are the bifurcate Killing horizons of maximally extended black 
hole geometries like the Kruskal extension of the Schwarzschild even including the non-extreme 
charged rotating case. 

For our purposes it is important to notice that, in the case of a bifurcate Killing horizon, 
any neighbourhood of a point on which has empty intersection with the bifurcation sur- 
face 23, satisfies the local general hypotheses stated in definition 12.11 It is very remarkable for 
physical applications and for our subsequent discussion in particular, that such a result can be 
partially reversed as established by Racz and Wald [RW921 IRW96) , Indeed, if the local general 
geometric hypotheses are fulfilled on !K, for a sufficiently small 0, the spacetime outside can 
be smoothly deformed preserving the geometry inside and extending K and "K to a whole 
bifurcate Killing horizon in the deformed spacetime. Thus, when studying local properties, the 
bifurcation surface 23 can be "added" also to those spacetimes without bifurcation surface, as 
is the case of black holes formed by stellar collapse. One can therefore take advantage of the 
various technical properties of the bifurcate Killing horizon as we shall do in the rest of the paper. 

2.3. Killing vector, geodesical distance in and coordinates V, U, x 3 , x 4 . Let us 

focus on a relevant coordinate patch |KW91] defined in a neighbourhood of "K+ for a bifurcate 
Killing horizon generated by a Killing vector field K (similar adapted coordinates exist for IK_ ) . 
Let U denote an affine parameter along the null geodesies forming "K + fixing the origin at the 
bifurcation surface 23. Each point p G is thus determined by a corresponding pair (U,s), 
where s G 23 denotes the point which is intersected by the null geodesic generator through p. 

We shall now extend those coordinate system to a neighbourhood of "K + . To this end, for 
each point q G 23, let us indicate by re the unique future-pointing null vector which is orthogonal 
to 23 and has inner product —1/2 with -^j. We extend re on all of "K + by parallel transport 
along the null generators of "K + . Let V denote the affine parameter along the null geodesies 
determined by re, with V= on 23. It is clear that (V, U, s) determines a point in a sufficiently 
small neighbourhood of Jf + . We are thus in a position to introduce the sought coordinate 
patch. If (x 3 ,x 4 ) denote coordinates defined on a open neighbourhood (in 23) of a point in 23, a 
coordinate patch (V, U, x 3 ,x 4 ), we call adapted to turns out to be defined in corresponding 
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Figure 1: The thick black line is part of the bifurcation surface. The region painted in light grey 
corresponds to the horizon !K, while the dark-grey region represents the open set 0x< 

open neighbourhoods of points on 5f+. In these coordinates: 

1 1 4 

g\ %+ = — dU ®dV dV ® dU + V] hij(x 3 , x 4 )dx l <g> dx ] , (1) 

2 2 . . 

where metric h is that induced by g on B. /i is thus positive definite and it does not depend on 
V, U. We stress once more that, in view of Racz-Wald's result, such a coordinate system always 
exits in 0, provided the local geometric hypotheses hold in 0, regardless the actual geometry of 
the spacetime outside 0. 

In the rest of the paper, referring to this geometric structure in 0, we shall employ the 
following notation. We shall indicate by §v,u the cross section of at V, U constant. Moreover, 
s(p) € S will be the point with coordinates (x 3 ,x 4 ) when p 6 has coordinates (V, U, x 3 , x 4 ) 
(similarly (V' , U' , x' 3 , x' 4 ) will denote the analogous set of coordinates for p' £ 0). The set 
Gs(p, V, U) C §v,u is defined as the set whose image s(Gs(p, V, U)) C S coincides with the open 
/i-geodesical ball centred on s(p) with radius S. 

We shall denote by o~(p,p') the squared ^-geodesic distance, taken with its Lorentzian sign, 
between any couple of points p,p' contained in some g-geodesically convex neighbourhood. Fi- 
nally, we denote by (t(s,s ! ) the squared /i-geodesic distance between points s,s' in some h- 
geodesically convex neighbourhood contained in S. 

Making use of the preceding definitions and the introduced notations, we are going to present 
the following useful proposition, which is also based on achievements in [K W91| . 

Proposition 2.1. Let C M be a set on which the local general geometric hypotheses hold, 
and let Oi be another open set such that Ox C is compact. Define a coordinate frame adapted 
to !K + on (as said at the beginning of section \2.3\) so that p £ has coordinates (V, U, x 3 ,x 4 ) 
with p £ "K iffV = 0. The following holds. 
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(a) In 0, the decomposition K = K 1 -^ + K 2 -^j + K 3 -£^ + K 4 -£^ is valid and, if p G Oi: 

K 1 (p) = -nV + V 2 Ri(p), K 2 ( p ) = kU + V 2 R 2 (p), K\p) = VRi(p) ,i = 3, 4 , (2) 

where R\, R2, Ri are bounded smooth functions defined on 0\. 

(b) If Oi is included in a g-geodesically convex neighbourhood and if p,p' G Oi D Oi, i/ien 
cr(p,p') = e(s(p),a(j/)). 

(c) // Oi is as m f&j, there exists a 5 > suc/z i/j.at ; /or every fixed p E 0±, the smooth map 
Gs(p, V, U') 3p'4 a(p,p') has vanishing gradient (with respect to the coordinates of s(p') ) 
in a unique point q(p,V ,U') attaining its minimum there. In particular s(q(p,V ,U')) = 
s(p) ifp G Oi. 

(d) For p and q = q(p, V , U') as in (c): 

a(p, q) = £(s(p), 8 (q)) - (U - U')(V - V) + R(p, V', V) , (3) 

where R(p, V , U') = AV 2 + BV' 2 + CVV' , for some bounded smooth functions A, B, C of 
p, V',U'. 

Proof, (a) V a K b + V b K a = and V a (K b K b ) = -2nK a on Oi imply that V K K = kK on Oi, so 
that K \m= kU j^j because U is an affine parameter and K vanishes on H> where U = 0. If x a 
is any of x 1 = V, x 2 = U, x 3 , x 4 , exploiting the parallel transport used to define the coordinates, 
we have 

Til tj£= T" 2 f 5{ = T22 \%= Tfx \-K= \-K= \ji= ^ll \%= ^la \%= . (4) 

Above, the fifth one is equivalent to g(-^j, V_s_ \-k= and it arises from g(-^j, gfa) = — n a 
and V_s_t^t = on Oi, the seventh one can be proved similarly. Next, taking the first-order 

dU OU 

Taylor expansion in V of both K 1 and g ab about V = 0, we have, for some smooth function 
V 2 Ri(p) bounded in view of the compactness of Oi, 

_ K^W + V (^«^« +*bc + V"*<p) . (5) 

From ^ = -5 M r^ - 5 ad r^ d and and ©, exploiting Q and if 1 = on 5C, the identity 
([5]) simplifies: The last derivative vanishes and K 1 = —V^r \ji +V 2 R 1 {p). Furthermore, 
V 2 Kx + ViK 2 = evaluated on Oi and using gj) leads to ^ fj { ;= [jf, so that: 

= +V 2 Ri(p) = V (g la \n +K a \x +V 2 R l {p) . 

The last derivative vanishes in view of 0), ([1]) and the know identity = gbdFia + 9adF d c b- 
Therefore the first identity in ([2]) holds in view of ([T]) and K\% R = ull-^j. The second one can 
be proved with the same procedure noticing that the Killing identity Vyifi = 0, on Oi, becomes 
-Mr \jc = i n view of flD). The last identity in ([2]) is obvious. 
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(b) Since the geodesically convex neighbourhoods form a base of the topology and the 
projection n : p i— > s(p) is continuous, if Oi is chosen to be sufficiently small, we have that Oi is 
contained in a geodesically convex neighbourhood while, at the same time, vr(Oi) is contained 
in a /i-geodesically convex neighbourhood in "B. Without loss of generality, we can further 
assume that the latter is included in a ^-geodesically convex neighbourhood of M. Thus cr(p,p'), 
a(s(p)), s(p')) and £(s(p), s(p')) are simultaneously well-defined if p,p' £ Oi PlIK for a sufficiently 
small Oi. We notice that, a(p,p') is invariant under the action of the Killing isometry. Hence, 
for any p,p' £ "K n Oi we get the identity c(p,p') = a(s(p), s(p')) taking the limit towards 23 of 
the flow generated by the Killing field K applied to p,p'. Finally a(s(p), s(p')) = £(s(p), s(p')) 
because 23 is totally geodesic as it can be proved by direct inspection. 

(c) Let (V, U,s)=p and (V , U' , s') = p' . Whenever both points p and pi are contained on the 
horizon, namely V = V' = 0, the thesis holds in view of (b) and the fact that £(s, s') is positive 
definite, with positive-definite Hessian matrix in the coordinates x /3 ,x /4 of s' . Furthermore, in 
this case s(q(p,0,U')) = s(p). By continuity, that Hessian matrix remains positive definite if 
p,p' stay close to JC, so that, any zero q(p, V , U') of the x' 3 , x /4 -gradient of §>v',U' 3p'H> a(p,p') 
determines a minimum of a(p,p'). Taking the Taylor expansion of V x iiO~(p,p') (i = 3,4) centred 
on a point in Jfx J{ with respect to all the coordinates of p and p', the equation for q(p, V' , U') 
can easily be handled by exploiting Banach's fix point theorem, proving the existence and 
the uniqueness of q(p,V ,U') for p £ Oi sufficiently shrunk around IK, and p' varying in a 
neighbourhood Gg(p, V, U') of (0, U' , s) in §v',U'- We recall that G$(p, V , V) is the pre-image 
through §v',U' 9p'^ s ip') of a geodesic ball on 23 centred on s(p). The compactness of Oi and 
a continuity argument assures that S > can be chosen uniformly in p. 

(d) Keeping U, U', x s ,x 4 fixed, the expansion ([3]) is nothing but the first-order (V, V)-Taylor 
expansion of a(p, q(p, V , U')) at V = V = 0, paying attention to the fact that the coordinates 
x' 3 , x' 4 of q(p, V', U') depend on V through p. □ 

3 Correlations across the Killing horizon 

3.1. General outlook. We wish to compute the correlation functions of a real scalar quantum 
field, for field observables localized in a region intersecting a Killing horizon. Thus, we assume 
that a quantum system is described by a corresponding *-algebra A generated by the unit I and 
field operators $(/) for all / £ C^°(M). For our purposes only minimal properties of such 
a quantum field theory are necessary, namely that the (abstract) field operator $ is smeared 
with compactly supported smooth functions / £ C^°(M) and furthermore that: (i) the map 
/ i y $(/) is linear, (ii) <&(/)* = $(/), and (iii) [«£(/), <&(<?)] = when supp(f) and supp(g) are 
causally separated. VFe stress that do not assume that $ satisfies any specific field equation, so 
that, in principle our approach may encompass interacting fields. 

If u is an (algebraic) state on A, the correlation functions we are interested in are the bilinear 
functionals that map real smooth functions /, /' to (/)<&(/')). We shall specify shortly the 
form of the test functions /, /'. For the moment we only say that their supports are taken in 
a open region containing a Killing field K and which satisfies the local general geometric 
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hypotheses. We finally assume that can be covered by coordinates adapted to "K + (defined at 
the beginning of section 12. 3p with U and V increasing towards the future. 

Later we shall restrict to a subregion Oi as in (b) of proposition 12.11 because we want to 
use the expression ([3]) for the geodesic distance. The region Oi when considered in coordinates 
(V, U, s), can always be taken of the form (—a, a) x (b, c) x S, where S C H> is an open relatively 
compact subset. Notice that, shrinking Oi around the region of 'K n Oi determined by (6, c) x S 
means taking a > smaller and smaller. 

In the sufficiently small neighbourhood Oi one finds g{K,K) = k 2 UV + 0(V 2 ) in view of 
(PQ), (a) in proposition 12.11 and dyg22\'K = (arising from r| 2 |j{ = 0). As a consequence, K 
turns out to be spacelike in O s = {p £ Oi | V(p) > 0} and timelike in t = {p S 0\ \ V(p) < 0}. 
Referring to stationary black holes, Oi can be interpreted as a sufficiently small region around 
a point on the future horizon, the only horizon existing when the black hole is produced by 
collapsing matter. There, O s is part of the internal region, containing the singularity, while Ot 
stays in the external region, stationary with respect to the Killing time associated to K. In this 
way, a notion of energy related to K can be defined, in Ot at least, and we will take advantage 
of it shortly. 

We recall that, in the GNS representation of a state u on the *-algebra of field observables, 
the expectation value of the product of two fields u) ($>(/) <&(/')) (with /, /' real) is equal to 
($>(f)fy UJ \&(f')fy UJ ) , where ^ is the cyclic vector and $(/) is the field operator represented as a 
proper operator on the GNS Hilbert space. Hence, up to normalization, \uj (<!>(/) $(/')) | 2 can be 
interpreted as a transition probability between the states $(/) x I' CJ and In this sense, 

when / and /' are localized on the opposite sides of the horizon, the regions S and Ot, the 
correlation function u (<I>(/)$(/ / )) provides a measure of the transition probability through the 
horizon. 

Inspired by the ideas proper of the scaling-limit procedure [HNS841 !Bu96l IBV95j . in order to 
obtain the leading order to that probability, we shall consider some sequences of smearing 
functions f\ and /{ whose support become closer and closer to the horizon "K in the limit 
A — > + . Thus we are going to compute the limit: 



where f\ and /( are smooth functions supported in Oi whose supports become closer and closer 
to the horizon as long as A — > + . 

In contrast to the discussion presented in |HNS84| . where the scaling limit towards a point is 
analysed, with the proposed procedure we are instead considering scaling limits towards a null 
surface, so that the results presented in |HNS84j cannot be automatically applied to the present 
case. 

Since only the short distance behaviour of the two-point function of the state is relevant for 
our computation, we select the class of allowed states looking at their ultraviolet features. We 
assume that the two-point function of uj is a distribution of D'(M x M) defined as 



lim . u (<&(/aWa)) 
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where the integral kernels ui e have the following form 



u} e (x, x') 




+ w e (x, x ) 



(6) 



whenever the test functions are supported in a fixed, relatively compact, geodesically convex 
neighbourhood. In the previous expression, a e (x,x') = a(x,x') + 2ie(T(x) — T(x')) + e 2 and T 
is a fixed (arbitrarily chosen) time function |KW91| . The smooth strictly-positive function A is 
the so-called Van Vleck-Morette determinant [Wa94^ IKW91| . We finally assume that w e has a 
"less singular" behaviour with respect to that of a t (x, x') -1 in the sense we are going to specify. 
An important case is: 

w e (x, x) = v(x, x) lncr e (x, x ) + w(x, x) for some fixed smooth functions v, w. (7) 

With this form of uu e , the right-hand side of ([6]) is a straightforward generalization of the short 
distance structure of the two-point function of Minkowski vacuum for a Klein-Gordon scalar 
field, if one also suppose that the field $ satisfies the Klein-Gordon field equation □<& + = 
(where V : M — > R is any fixed smooth function). In particular, ©-(JT]) are fulfilled by all the 
quasifree states of Hadamard type [Wa94] that are defined by ([6])- ((7|) with a further requirement 
on the form of v related to the Klein-Gordon equation. Those states are supposed to be the 
most significant states for free QFT in curved spacetime [W&94] and are very often employed in 
the rigorous description of thermal properties of quantum fields in the presence of black holes 
[KW91, FH90, DMP09]. However, since we do not need to consider free fields or any precise 
field equation, we further relax the requirements, adopting ([6]) but, in place of ([?]), assuming a 
weaker pair of requirements (in the following p = (V, U, s), p' = (V , U', s')): 

(Wl) w e (p,p') — > w'(p,p') as e — > + , almost everywhere in (p,p'), for some function w' and 
w e is e-uniformly bounded by a locally M 2 -integrable function; 

(W2) w'(V, U, s, V, U', sf) -> w"(U, s, U', s') almost everywhere in (U, s, U', s') when (V, V) -»• 
(0, 0) for some function w" on J{ 2 , and w 1 is (V, V')-uniformly bounded by a locally 5{ 2 -integrable 



Notice that (Wl) and (W2) are satisfied by w t whenever it satisfies (|7|) and thus they are 
valid for Hadamard states in particular. 

To conclude, few words about the precise construction of the functions f\, f[ are necessary. 
Let / and /' be some smooth functions with compact support contained respectively in the 
regions O s and Of The associated functions fx, fx are defined as 



f x (V, U, x 3 ,x 4 ) := -f -, U, x 3 , x 4 , f x (V, U, x 3 , x 4 ) := -f -, U, x 3 , x 4 , A > . (8) 



As in other scaling-limit procedures, the pre-factor A is introduced in order to keep the result 
finite. Finally, to avoid divergences due to zero- modes in the limit A — > + , we assume that /, /' 



function. 



1 




1 
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are of the forrrjl: 

dF OF' 

f=Qy> f = W forfixed ^ f ' eC o°(Oi)- (9) 

Alternatively, sticking to generic smooth compactly supported functions / and /', the divergent 
contribution of zero- modes has to be subtracted at the end of the computations. With the as- 
sumptions (|HJ) and ©, the A — > + limit of to (<£(/>)<J>(/^)) is precisely the notion of scaling limit 
of lo (dv&(x)dv&(y)) towards the horizon we shall employ. The result of such limit represents 
the first contribution to the sought transition probability in an ideal asymptotic expansion for 
small A. 



3.2. Computation. We are in a position to present the most important result of this paper. 
The proof of the following theorem exploits techniques similar to those employed in the appendix 
B of [KW91]. However, our result differs from those presented there because we are not interested 
in studying the restriction of the states to the horizon in the symplectic approach, while we intend 
to compute the scaling limit of the state in the AD smeared formalism. The relevant part for 
our computation arises from the (V — V) contribution to a in ([3]) whereas, in [K W91| . only 
the (U — V) contribution plays a relevant role. In this respect, there are some similarities to 
the analysis performed in [FH90], although here neither we restrict ourself to the spherically 
symmetric case nor we consider any equation of motion for the quantum field, nor we suppose 
that the Killing structure extends up to the future null infinity, as done in [FH90 . 

Theorem 3.1. Assume that the general local geometric hypotheses as in Def \2.1\ hold in 
(covered by coordinates adapted to "K+), suppose that Oi C is a sufficiently small open neigh- 
bourhood of a point on "K with 0\ C compact. Assume also that the state u) has two-point 
function given by a distribution as in ([6]) that verifies the requirements (Wl) and (W2) above 
(in particular, w t may have the form ^ and may be a Hadamard state for Klein- Gordon fields). 
If f,f are taken as in |2)j-(0J) and \x is the measure associated to the 2-metric on the bifurcation 
surface H, it holds 

lim u (HfxMf'x)) = lim - JL f HV ;^ ^'y':^ S) dUdVdU'dV'd^s) . (10) 



A 



Proof. We start by focusing on the contribution of the most singular part of the two-point 
function ([6]), that is the iterated limit: 

£:=limlim / ^b,^)/ > M&/) 

where dp is a short-cut for the measure induced by the metric. Fixing 5 > and p, V',U', consider 
a neighbourhood G$(p, V , U') C §v',U' as m ( c ) of proposition 12.11 and define a smooth map 23 B 



2 Indeed, given / 6 C§°(0i), an F G C£°(0i) with / = ff exists if and only if / R f(V, U, x 3 , x 4 )dV = on Oi, 
namely, if and only if f(-,U,x 3 ,x 4 ) has no zero modes referring to the V-Fourier transform. 
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s h-» xs( s , V, U',p) > with support completely included in Gg(p, V', U') and x&i s ' V 1 U' ,p) = 
1 for < yj X(s(p), s') < | + ^y/\(s(p), s(q(p, V, U'))). In view of the smoothness of all 
considered functions it is possible to arrange these functions in order that (s' ,V' ,U' ,p) —> 
X5(s' ,V ,U' ,p) is jointly smooth. Finally we can decompose the integral in (fTT]) as: 

A 1/2 (p,p')f\(p)f' x (p') , , , , , / , , iW JXflJTTtA '\ 

A 2 ( 7T Xs{s,V,U,p) a/ \ det g{p')\dV dU ds 

4Tr^a e {p,p') J 

^MhWf'xW) (1 _ y, ^ }) V| det (12) 

Oi \JQl ^ 2 a t (p,p') I 




+ 



where ds' = dx' 3 dx' 4 . Let us start from the second integral. As consequence of the compactness 
of Oi, the continuity of a and (b) in proposition 12.11 we have that, for a fixed r\ > and for 
p,p' in a sufficiently small 0±, y cr(p,p') > rj/2 if ^/£(s(p), s(p')) > r\ . We stress that the limit 
in A — > + in ([IT]) allows us to take Oi as small as we need. Thus, by definition of x<5 an d of 
Gs(p, V, U'), the denominator a t= o(p,p') > 5 2 /16 when 1 — xs ^ 0. Hence, the integrand of the 
second integral in (|12p is jointly smooth in all variables including e, even for e = 0. Then, in view 
of Lebesgue's dominated convergence theorem, the limit in e can be computed simply taking 
e = in the integrand. For the same reason, after changing the integration variables (V, V) to 
(XV, XV'), the subsequent limit as A — > + can be computed under the sign of integration. The 
resulting integral vanish because, in the (V, V) variables, it is noting but the integral of the V, V 
derivative of some compactly supported smooth function. Thus only the former integral in (|12p 
may survive the limits in (jlip . Let us focus on that integral. Making use of (c) in proposition 
12.11 in each set G$(p, V , U') C §>v',U' we define the function p(p') > such that: 

o{p,p) = p(p') 2 + cr(p, q{p, V', U')) . 

In view of (c) in proposition 12.11 the pair p, 6, where 6 £ (—ir,ir) is the standard polar angle in 
geodesic polar coordinates centred on q(p,V ,U'), determines an allowable local chart for any 
p' G Gs(p, V , U') (see also the appendix B of |KW91j ), that is smooth barring the usual conical 
singularity for p = 0. Notice that, due to the last statement in (c) of proposition 12. 1\ when 
p € !K and V' = 0, p coincides with the standard geodesic radial coordinate centred on s(p) G !B. 
In the following we shall employ that coordinate system in each Gg(p,V',U'). Making finally 
use of (d) in proposition 12.11 and choosing T = (U + V)/2, we can re-arrange the former integral 
in CG]) so that: 

T = v v f &' 1/2 (p,p')fx(p)f'x(p') <¥dp 

ATo e "o+i p 2 -(V-V'-ie)(U-U'-ie) + R(p,V',U') 4vr 2 ' [ ' 

where, A n / 2 (p,p') := A 1 / 2 (p, p')x&{p, p') ■ From now on we shall assume that the integral in p' 
is performed before that in p. Using this coordinate system the integral in the right-hand side 
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of (p73|) can be rewritterU as 



jA n /\p,p')f x (p)f' x (p')-?-ln {p 2 + a(p, q{p, V, U'))) dpdOdU' dV dp 

where det g is the determinant of the metric in the coordinates p, 9, V, U' , which parametrically 
depends on p. Notice that the domain of integration in p is bounded by the support of the 
function Xs{PiP') embodied in A'. For V = V = 0, the metric takes the form ([1]) on ¥> which 
does not depend on U, U', V, V any more while R vanishes. By direct inspection one sees that 

^ dct ^ is continuous, tends to 1/2 whenp — > p' in "K and its p-derivative is continuous for p ^ 
it being however bounded there. If A' A , R\, detgA> dp\ are respectively defined as A', R, detg 
and dp with V and V' rescaled by A, changing coordinates (V, V) — > (XV, XV) the integral in 
the right-hand side of (|13p can be rearranged as 



Jd v F(p)d v ,F'(p')A't /2 (p,p')-?- ln[p 2 - (XV - XV - ie)(U - V - ie) + R x (p, V, U')] 

^^^ dpdOdU'dVdpx . 

Leaving unchanged the remaining integrations, we can first integrate by parts in the polar 
coordinate p. We obtain two boundary terms, which are integrals in the remaining variables 
evaluated respectively at p = and p = po sufficiently large, and an integral in all the variable 
including p. When taking the limits as e — > + and A — > + , concerning the integral representing 
the boundary term at p = po, we can pass both limits under the sign of the integration by 
straightforward application of Lebesgue's dominated convergence theorem. The result is that, 
in the integrand, only dyF and dy>F' depend on V, V , hence, performing the integrations in V 
and V both integrals vanish because F and F' are of compact support. The remaining boundary 
term leads to the limit: 

L= lira lim -2vr [ A' 1/2 (XV,U, s, XV ,U' , s)d v F(V,U, s)d v/ F' (V' ,U' , s) 

A^0+ £->0+ J 



8tt 2 P 



dU'dVdpx 

p=0 



Above, the factor 2ir arises by the integration in 6 at p = 0, and we have safely replaced e with 
Ae in the integral in view of the order of the limits. Notice that, thanks to (d) of proposition 
12.11 \X^ 1 R\(p, V , U')\ < CX for some constant C, when p,p' £ Oi and A G [0, Ao). Furthermore, 
from the (pQ) of the metric on "K we get that dp x = |(1 + XV z)dU dV dp{s) for some smooth 
function z = z(V,U,s), where dp is the measure associated to the 2-metric h on S. Thus the 
term proportional to log (A) can be dropped as it gives no contribution to the final result because 
go, Ao do not depend on V and the apparently divergent term results once more in the integral 



3 We henceforth assume to cut the complex plane along the negative real axis to define the function lnz. 
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of a derivative of a compactly supported smooth function. The limits of the remaining terms 
can then be computed, in the given order, exploiting Lebesgue's theorem. We eventually obtain: 

/"o r,o rv/ i it, „,„„ rH,\ dUdVdU'dV'duis) . A , 
L = - J dvFdv>F'(iTTXE + XA + -iKXE + XA-+ln\V-V'\\U-U'\) — (14) 

where E± is the subset of Oi with, respectively, (V — V')(U — U') ^ and A± is the analogue 
with, respectively, U — U' ^ 0, and xs is the characteristic function of the set S. We have also 
used the fact that dp\ becomes ^dlldVdp(s) for A = in view of the form ([T]) of the metric 
on "K and that A' = A = 1 when s = s' and V = V = (it follows from A(p,p) = 1 and, 
since A is invariant under isometries, using an argument similar to that employed to prove (b) 

of proposition 12. ip . For the same reason, in view of the meaning p, ^/ dct ^ 9x ^ _ >. 1/2 for p — > 



when A = 0, when working in coordinates p, 9, V, U. The integral in (|14|) can equivalently be 
re-written introducing another e-prescription as: 

L= lim — f d v Fd v/ F' ln(-(V -V' -ie)(U -U'))dUdVdU'dV'dp(s) . 

e^0+ loir J 



Eventually, integrating by parts, we obtain (|10p concluding the proof, provided the contribution 
of the term w t in ([6]) to the left-hand side of (|10h vanishes. Indeed this is the case. Taking 
the requirements (Wl) and (W2) into account, applying Lebesgue's dominated convergence 
theorem, changing variables V —> XV , V — > XV and exploiting Lebesgue's theorem again, we 
have: 

f 11 1 1 f 11 1 1 dF dF' , , 

lim lim / f\(p)f\(p)w t (p,p) dpdp = / w (U,s,U',s') — ——dUdVdp(s)dU'dV'dp(s') = 
A->0+e^0+J J aV oV 

in view of the fact that the result of the integrations in V and V vanishes because F and F' 
are smooth with compact support. □ 

3.3. The correlation functions and their thermal spectrum. As is known (e.g., see 
[Wa84]) a timelike Killing vector field on the one hand provides a natural notion of time, which 
is nothing but the parameter of the integral lines of the field. On the other hand it gives a 
natural notion of conserved energy for fields and matter propagating in the region where the 
Killing vector is present. We are interested in computing the energy spectrum of the correlation 
functions uj(<fr(f\)<&(f' x )) seen by an observer that moves along the curves generated by the 
Killing field K and computed with respect to the associated Killing time. More precisely, 
exploiting Theorem 13. 1\ we intend to compute that energy spectrum in the limit of test functions 
squeezed on the local Killing horizon. As the supports of the test functions are infinitesimally 
close to the horizon, we have to focus on what happens for V ~ 0. Therefore we truncate every 
component of the formula ([2]) for the Killing vector at the dominant order in powers of V and 
we make use of the right-hand side of (I10p as definition of correlation two-point function. Now 
r is the Killing time, namely the integral parameter of the curves tangent to K. In the said 
approximation, the first identity in the right-hand side of j2j) implies: 

V{t) = -e~ Kr for V < (that is in t ) and V(t) = e" KT for V > (that is in S ), (15) 
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up to an additive constant in the definition of r which in principle could depend on the integral 
curve. Our choice is coherent with the standard definitions of r in Schwarzschild spacetime 
where r is the Killing time in the external region. Indeed we recover those cases in the limit 
where the r-constant 2-surfaces are close to the Killing horizon. We now examine two situations. 

(a) Both the supports of f\ and f' x stay in Of 

In that case, thinking of the functions F, F' as functions of r, r' instead of V, V in view of (|15p . 
we can re-arrange the found expression for the correlation function as 

i^WAWM = mg. / ^g^^_ iTMiu ^ {x) , {16) 

where we have used the fact that the functions F and F' are compactly supported by construction 
even adopting the new coordinate frame. It is known that, in the sense of the Fourier transform 
of the distributions, f R ^ ( sin h(f )+»()+ ) 2 = e ^-e-^ E ( e ' g '' see the a PP endix of |DMP09j). 

That identity and the convolution theorem lead to 

where fin = 2tt/ k is the inverse Hawking temperature and we have defined the smooth function: 

F(E,U,x) := [ -^=e- iET F(T,U,x), 

that, for E — > ±00, vanishes faster than E~ n , n = 0, 1, 2, . . . uniformly in the remaining coor- 
dinates. The thermal content of the found correlation function is manifest in view of the Bose 
factor (1 — e~ /3llE )~ 1 where the Hawking temperature 1/(3h takes place - fixed with respect to 
a (generally arbitrary) choice of the scale necessary to define the Killing field K (see the remark 
below). 

(b) The support of fx stays in O s , while that of f' x stays in O t ; Tunnelling processes. 

As previously remarked, up to normalization, \uj(<&(f\)$>(f x ))\ 2 can be interpreted as a tunnelling 
probability through the horizon. Employing (|15p once more, we end up with: 

iimk>($(/ A )*(/j)) = £ J w^w **,**?***) ■ m 

As expected from the fact that, in this case, the support of f\ is always disjoint from the support 
of f' x , we can directly pass the limit e — > + under the sign of integration simply dropping ie in 
the denominator. Taking advantage of the convolution theorem, the final result reads: 

h W (*(/A)*«D) = Tcl 1 r ^^(r'^^^ ^z) dUdU 'M*) ■ (19) 
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Notice that, if the arbitrary additive constant defining r in t were different from that in O s , 
then a further exponential exp (icE) would take place in the numerator for some real constant 
c. In any case, the energy spectrum does not agree with the Bose law. However, considering 
packets concentrated around to a high value of the energy Eq, (I19p leads to the estimate for the 
tunnelling probability: 

Jim \uj($(fxMf x ))\ 2 ~ const. E 2 e~^ E ° , 

A— s>0 

in agreement with the ideas in [PWOOl £NVZ05j . It is nevertheless worth remarking that the 
interpretation of E as an energy is questionable for the packet in the internal region S since 
the Killing vector K is spacelike therein. 

Remark@. It is worth stressing that if K is a (future-directed) Killing field, for every constant 
c > 0, cK is a (future-directed) Killing field again. In general, there is no way to fix the scale 
c. This arbitrariness in the definition of K enters the definition X7 a (Ki,K b ) = —IkK" 1 of the 
associated surface gravity k and, in turn, it affects the definition of the Hawking temperature 
Tjj = k/(2tt). In the case of an asymptotically flat (say, Schwarzschild) black hole, one uses the 
asymptotic behaviour of the Killing isometries to normalize the "Schwarzschild time" parameter 
to correspond to "ordinary time" at infinity (and of course energy is also defined with respect 
to infinity), fixing the scale c. However, if one is assuming only the existence of a local horizon 
and not necessarily any connection with any "infinity" region, in the absence of further physical 
requirement, there is nothing to set the scale. So, while it may be correct to say that one has 
a thermal spectrum in the found expressions for lim>_>o w (3K/a)3K/a))> it does not seem very 
meaningful to claim that this corresponds to any particular value for the Hawking temperature 
in the general case. 



4 Conclusions 

In this paper we have analysed the correlation functions and the tunnelling amplitude through a 
Killing horizon for quantum states in a certain, physically relevant, class of states that includes 
the states of Hadamard type. Although the equation of motion does not play a relevant role, the 
results are in particular valid for scalar Klein-Gordon particles. The computation is performed 
in the limit of test functions squeezed on the Killing horizon. The considered local Killing 
horizon with positive constant surface gravity may be a part of the complete horizon of a 
black hole, including non-static black holes as the non-extremal charged rotating one, or it 
may just temporarily exist in a finite region. The considered states are generally not required 
to be invariant with respect to the isometry group generated by the Killing field. We have 
established that, in the limit of test functions sharply localized on the opposite sides of the 
horizon, the correlation functions have a thermal nature, namely they have a spectrum which 
decays exponentially as exp{— /Snawking^} for high energies. The energy E and the inverse 

4 The authors are grateful to one of the referees for this important remark. 
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Hawking temperature are denned with respect to the Killing field generating the horizon (see 
however the remark at the end of the previous section). This achievement is in agreement with 
the result obtained in other recent papers, although here it is obtained in the framework of 
the rigorous formulation of quantum field theory on curved spacetime. While, in the mentioned 
literature quantum-mechanical approaches are exploited, leaving unresolved the issue concerning 
the strongly ambiguous notion of quantum particle in curved spacetime. 

Furthermore, we have also established that, when both test functions are localized in the 
external side of the horizon, a full Bose spectrum at the Hawking temperature arises in the 
expression of the correlation functions. In both cases the computation is completely local, i.e. 
the nature of the geometry at infinity does not matter and the results do not depend on the 
employed states provided they belong to the mentioned class. These results give a strong support 
to the idea that the Hawking radiation, that it is usually presented as a radiation detected at 
future (lightlike) infinity and needs the global structure of a black hole Killing horizon, can 
also be described as a local phenomenon for local geometric structures (local Killing horizons) 
existing just "for a while". A fundamental ingredient in our computation is the fact that the 
nonvanishing surface gravity is constant on the Killing horizon. Thus, we could exploit the result 
m [RW92] and, in turn, some technical constructions of [KW91]. Even if this requirement can 
easily be physically interpreted as the geometrical description of the thermodynamic equilibrium, 
it would be interesting to consider, from our viewpoint, the case of a black hole in formation, 
where there are no Killing horizons at all. The latter situation has already been investigated, 
at least in the presence of spherical symmetry, as in [D NVZZ07t IHDV NZ09], In those papers 
the WKB approach as well as the theory of Kodama-Hayward and the associated notion of 
dynamical horizon are exploited. As a preliminary comment, we notice that the computation of 
the scaling limit towards the horizon does not seem to require the presence of a proper Killing 
horizon, which could be replaced by some more generic null hypersurface. 
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